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Abstract 

The  technique  of  Picosecond  Time-Resolved  Laser-Induced  Fluorescence  (PITLIF)  was 
expanded  to  enable  two-point  hydroxyl  (OH)  concentration  time  series  measurements  of 
turbulent  jet  flames.  Measurements  were  obtained  with  two-point  separation  along  the  radial 
direction  of  the  jet  (Ar),  the  propagation  direction  of  the  horizontal  laser  beam,  and  axial  two- 
point  separations  along  the  jet  axis.  As  the  jet  configuration  was  statistically  axisymmetric, 
statistical  analyses  along  the  azimuthal  direction  could  be  considered  insignificant,  allowing 
complete  coverage  of  the  jet  flame.  Investigation  of  axial  two-point  statistics  required  a  different 
laser  configuration  than  that  of  the  previous  radial  two-point  system:  either  a  single  vertical  laser 
beam  or  two  horizontal  beams  with  adjustable  vertical  separation.  The  combustor  used  in  this 
study  had  a  closed  base  that  did  not  permit  vertical  alignment  of  a  laser;  hence,  the  two-beam 
configuration  was  adopted.  With  the  installation  of  our  new  diode-pumped  laser  system,  which 
generates  nearly  ten  times  more  UV  power  than  before,  the  laser  could  be  split  into  two  beams 
while  still  providing  a  sufficient  signal  level  at  each  probe  volume. 
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Technical  Discussion 


Picosecond  Time-Resolved  Laser-Induced  Fluorescence  (PITLIF) 

Quantitative  measurements  of  radical  species  concentrations  in  flames  are  required  for  an 
understanding  of  important  interactions  between  fluid  mixing  and  chemical  reactions  [1]. 
Advances  in  laser-based  techniques  have  made  non-perturbing  measurements  of  such  species 
concentrations  possible  in  recent  years.  In  particular,  laser-induced  fluorescence  (LIF)  possesses 
both  the  spatial  and  temporal  resolution  necessary  for  monitoring  radical  concentrations  in 
reacting  flows  (Barlow  and  Carter,  1994)  [2].  Unfortunately,  LIF  usually  suffers  from  an  inverse 
dependence  on  the  local  quenching  rate  coefficient.  Although  laser-saturated  fluorescence  (LSF) 
permits  measurements  of  species  concentrations  without  recourse  to  quenching  calculations,  the 
LSF  technique  requires  large  input  energies  that  demand  a  laser  repetition  rate  on  the  order  of  10 
Hz.  This  repetition  rate  can  provide  mean  concentrations  and  probability  density  functions 
(PDFs),  but  the  10  Hz  repetition  rate  is  clearly  too  slow  to  resolve  species  fluctuations  in 
turbulent  flames.  High  repetition  rate  lasers,  on  the  other  hand,  can  be  utilized  for  measuring 
turbulent  species  fluctuations  but  the  quenching  rate  coefficient  must  be  obtained  within 
sampling  times  on  the  order  of  the  Kolmogorov  time  scale. 

As  a  consequence  of  the  above  limitations  for  the  usual  LIF  methods,  the  LIF  protocol 
chosen  for  this  research  program  utilized  a  mode-locked  laser  that  had  the  temporal  resolution 
and  irradiance  necessary  to  provide  direct  time-series  measurements,  thus  yielding 
autocorrelation  functions,  as  well  as  more  traditional  PDFs.  The  resulting  optical  technique, 
picosecond  time-resolved  laser-induced  fluorescence  (PITLIF),  operated  by  monitoring  the  LIF 
signal  using  a  novel  four-bin  integrated  photon-counting  procedure.  Three  sequential  bins  were 
used  to  resolve  temporal  decays  rapidly  to  obtain  fluorescence  lifetimes  at  frequencies  exceeding 
10  kHz.  The  remaining  bin  was  integrated  over  the  entire  fluorescence  lifetime  to  obtain  the 
integrated  fluorescence  signal.  The  quenching  data  from  the  first  three  bins  ultimately  could  be 
used  to  correct  the  integrated  data  from  the  remaining  bin,  thus  producing  a  temporal  record  of 
number  density  fluctuations  that  were  independent  of  the  local  quenching  environment.  From 
this  temporal  record,  quantitative  autocorrelation  functions  and  power  spectral  densities  (PSD) 
could  be  determined  corresponding  to  local  concentrations  in  the  flame.  The  resulting  integral 
time  scales  then  could  be  used  for  validation  of  theoretical  turbulence  and  combustion  models. 

Over  the  past  three  years  we  implemented  PITLIF  for  measurements  of  the  OH  integral 
time  scale  using  a  wide  range  of  turbulent  jet  flames.  From  such  data,  suitable  statistical 
information  was.  compiled  for  comparisons  to  predictions  from  large-eddy  simulations  (LES). 
In  addition,  we  developed  a  new  version  of  PITLIF  that  permitted  simultaneous  time  series  at 
two  nearby  points  in  turbulent  combustors.  Before  the  two-point  systems  and  results  are 
presented  in  this  report  a  review  in  some  detail  of  the  theory  and  basic  photon-counting 
instrumentation  for  PITLIF  is  described. 
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Theory  for  PITLIF 

The  basic  PITLIF  procedure,  which  produces  quenching-independent  measurements  of 
concentration,  can  be  understood  by  applying  the  radiative  rate  equations  to  a  simple  two-level 
model.  We  begin  with  a  laser  pulse  that  excites  molecules  from  the  ground  level  to  an  excited 
electronic  level.  All  of  the  initial  population  is  assumed  to  be  in  the  ground  level  before 
excitation;  moreover,  the  width  of  the  laser  pulse.  At,  is  short  compared  to  the  quenching  time. 
For  sufficiently  low  laser  power,  the  number  density  in  the  excited  level  immediately  following 
the  laser  pulse,  N%  (cm'3),  is  given  by 

N;  =  N'xWnto,  (1) 

where  N°  is  the  initial  number  density  of  the  species  of  interest  and  IVI2  is  the  average  rate 
coefficient  for  stimulated  absorption  (s'1),  as  computed  from  an  integration  over  the  exact  laser 
pulse  shape.  Subsequent  to  the  laser  pulse,  the  population  of  the  excited  level  decays  back  to  the 
ground  level.  In  the  absence  of  laser  irradiation,  absorption  and  stimulated  emission  can  be 
neglected  so  that  during  this  decay  the  number  density  of  the  excited  level  is  given  by 

^2(/)  =  ^i0/Ti2A/exp(-//rn)>  (2) 

where  rn  is  the  excited-state  (or  fluorescence)  lifetime  (s).  The  excited-state  lifetime  is  given  in 
terms  of  the  rate  coefficient  for  spontaneous  emission,  A2I,  and  that  for  quenching,  Q2],  by 

r»=K,  + £,)''•  (3) 

The  temporal  fluorescence  signal  depends  linearly  on  both  the  number  density  in  the  excited 
level  and  the  collection  volume,  Vc  (cm3).  Hence,  the  temporally  resolved  fluorescence  decay  is 

given  by 

S(/)=lP12A//l21FA/;exp(-//ra),  (4) 

where  S(t)  represents  the  output  of  the  PITLIF  instrument  in  photons/second.  While  temporal 
analysis  of  S(t)  yields  the  fluorescence  lifetime,  integration  gives 

OC 

Sy  =  \s(t)dt  =  iVuAtAllVTnN°l,  (5) 

0 

which  corresponds  to  the  integrated  fluorescence  signal.  All  of  the  terms  in  Eq.  (5)  can  be 
determined  through  calibration  except  rfl  and  N°.  For  our  two-level  system,  the  excited-state 
lifetime  can  be  determined  from  the  decay  represented  by  Eq.  (4).  Equation  (5)  then  can  be  used 
to  determine  the  initial  number  density  corrected  for  electronic  quenching.  For  a  multilevel 
system,  the  total  number  density  can  be  obtained  from  N°  by  using  the  appropriate  Boltzmann 
fraction.  In  either  case,  the  integrated  fluorescence  signal  is  proportional  to  the  concentration  of 
target  molecules  in  the  flame. 


Previous  PITLIF  Measurements 

We  previously  demonstrated  that  our  photon-counting  procedures  for  PITLIF  can  be  used 
to  monitor  fluctuations  in  OH  concentration  on  the  time  scale  of  turbulence  [3-6].  During  the 
last  three  years,  we  applied  PITLIF  to  measurements  of  OH  statistics  in  turbulent  non-premixed 
and  partially  premixed  flames  [7,  8].  While  making  such  measurements,  we  examined  PITLIF 
statistics  and  LES  simulations  for  OH  in  turbulent  non-premixed  flames.  More  recently,  we 
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developed  a  two-point  version  of  PITLIF  that  now  permits,  for  the  first  time,  determination  of 
spatial,  as  well  as  temporal,  statistics  in  turbulent  flames  [9,  10].  Further  details  of  the 
experimental  arrangements  and  salient  results  are  reviewed  in  this  and  the  following  sections  of 
this  report. 


Experimental  Setup  and  Procedure 

A  diagram  of  the  laser  system,  including  the  burner,  is  shown  in  Fig.  1.  A  Spectra 
Physics  Tsunami,  regeneratively  mode-locked,  Ti:Sapphire  laser  is  pumped  by  a  20-W,  Spectra 


Gat'd  dliortmlnatora 


Physics  argon-ion,  multi-mode  laser  and 
provides  a  pulse  repetition  rate  of  80 
MHz.  For  OH  measurements,  the 
resulting  IR  beam  has  a  temporal  pulse 
width  of  2.0  ps  and  is  frequency-tripled 
in  a  CSK  SuperTripler.  The  output 
beam  is  re-collimated  by  two  UV  lenses 
and  focused  by  a  20-cm  focal  length, 
5.0-cm  diameter  UV  lens  to  form  the 
probe  volume  above  the  burner 
assembly.  For  the  tripled  beam,  the 
spatial  resolution  is  200x200x40  pm1 
based  on  the  beam  diameter  (e'2)  and  the 
monochromator  entrance  slit  width  for 
the  detection  system.  The  laser  power  is 
24-40  mW,  which  results  in  an  average 
probe  volume  irradiance  of  ~2.5xl07 
mW/cm2.  The  Ri(ll)  transition  of  the 
(0,0)  vibronic  band  (306.5  run)  of  OH  was  chosen  for  excitation.  This  line  displays  an 
approximately  ±5%  Boltzmann  fraction  variation  over  the  relevant  temperature  range  of  1500- 
2500  K. 

The  OH  fluorescence  is  collected  at  a  90-degree  angle  from  the  incident  laser  beam  by 
using  two  14.1 -cm  focal-length,  5.1 -cm  diameter  UV  lenses.  The  wavelength  of  the  measured 
fluorescence  is  selected  by  use  of  a  0.25-m  monochromator.  An  adjustable  slit  at  the  entrance  to 
the  monochromator  limits  the  probe  volume  along  the  beam  path.  For  most  measurements,  the 


iw 

0  26  m  monoohromitor 


Figure  I.  Experimental  setup  for  single-point  PITLIF 


Figure  2.  Graphical  representation  of  the  PITLIF  method;  D2,  D3,  and  D4  are  each  3.5  ns  wide.  The  background,  B , 
is  typical  of  flame  emission,  and  the  amplitude,  A ,  is  proportional  to  concentration. 

Figure  3.  Wire  schematic  for  the  PITLIF  system.  N:  NIM  input  from  the  discriminator;  I:  start  signal  input;  O:  start 
signal  output;  Cl:  channel  advance  input;  CO:  channel  advance  output.  Unused  signals  are  50-0 
terminated,  and  10.0,  6.5,  and  3.0  represent  the  temporal  delays  (ns)  placed  in  the  respective  lines. 
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spectral  window  was  adjusted  to  give  a  total  bandwidth  of  10  nm  centered  at  309  nm.  A 
Hamamatsu  HS5321  photomultiplier  tube  (PMT)  detected  the  fluorescence  at  the  exit  plane  of 
the  monochromator.  This  PMT  had  a  risetime  of  700  ps  and  a  transit  time  spread  of  160  ps.  By 
biasing  at  -2450  V,  the  single-photon  pulse  height  was  enhanced  for  subsequent  leading-edge 
discrimination.  Figure  2  shows  representative  fluorescence  decays  with  fluorescence 
lifetimes  (r„)  of  1.5  and  4.0  ns,  along  with  three  areas,  D2,  D3,  and  D4,  specifying  three  of  the 
four  measurements  for  the  gated  photon  counting  system.  The  remaining  measurement,  Di, 
which  is  not  shown,  is  not  gated  and  measures  all  photons;  thus,  Di  represents  the  total  integrated 
fluorescence  under  either  curve  in  Fig.  2. 

A  detailed  wire  schematic  for  the  gated,  photon-counting  system  may  be  seem  in  Fig.  3. 
Briefly,  this  system  consisted  of  two  LeCroy  Model  4608C,  eight-channel  discriminators 
connected  to  four  EG&G  Ortec  logic-pulse  counting  boards.  Several  channels  on  the  first 
discriminator  were  used  to  convert  each  PMT  pulse  into  four  logic  (NIM)  pulses.  The  second 
discriminator  gated  three  of  the  four  NIM  pulses  (D2-D4)  from  the  first  discriminator  using  a 
photodiode  output  from  the  Ti: Sapphire  cavity.  Separate  cable  delays  for  each  of  the  three  PMT 
inputs  to  the  second  discriminator  ensured  that  each  gate  passed  the  appropriate  portion  of  the 
fluorescence  decay.  The  ungated  PMT-NIM  pulse  from  the  first  discriminator  and  the  three 
outputs  from  the  gated  discriminator  are  counted  separately  by  the  four  EG&G  pulse  counters 
shown  schematically  in  Fig.  3.  These  comprised  the  measured  bins  D2-D4  of  Fig.  2  and  the 
integrated  measurement  Dj.  The  four  photon-counting  boards  each  had  an  8192-channel 
memory  and  could  be  sampled  simultaneously.  Each  channel  acquired  counts  over  many 
thousands  of  laser  pulses  as  set  in  software.  The  maximum  sampling  rate  (channel  advance  rate) 
for  the  boards  is  500  kHz,  well  above  that  needed  for  turbulence  studies. 

The  fluorescence  lifetime,  r„,  and  the  flame-emission  background,  B,  can  be  found  from  the 
gated  counts  using  the  expressions  (Pack  et  al.,  1998)  [11]. 


ln[(D2-D3)/(DrD)] 


A t(C2- 1)  ’ 


(7) 


where  At  was  the  width  of  each  integrated  bin  (3.5  ns)  and  C  =  exp(-A//rfl ).  Following  Eq.  (4), 
the  concentration  of  the  excited  state  immediately  following  the  laser  pulse  was  proportional  to 
the  amplitude  of  the  exponential  decay,  A,  shown  in  Fig.  2,  and  could  be  determined  from 


A  = - 2 - 1 - ,  (8) 

exp(-r0/rflXl-C2Xl-C)rfl 

where  tQ  was  the  delay  between  the  laser  pulse  and  the  start  of  bin  D2.  This  amplitude  was  also 
directly  proportional  to  the  concentration  of  the  ground  state,  following  Eq.  (1),  which  in  turn 
was  directly  proportional  to  the  total  number  density  of  the  molecule  of  interest  multiplied  by  a 
temperature-dependent  Boltzmann  fraction.  Since  the  excitation  wavelength  could  be  chosen  to 
minimize  experimentally  the  variation  of  the  Boltzmann  fraction  with  temperature,  A  in  Eq.  (8) 
could  be  taken  as  proportional  to  the  total  number  density. 
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Two-Point  PITLIF 


Radial  Measurement  System  and  Results 

This  two-point  technique  is  capable  of  simultaneously  monitoring  OH  time  series  at  two 
locations,  hence  allowing  the  characterization  of  spatial  structures  in  turbulent  flames  in  addition 
to  their  temporal  behavior.  Several  refinements  have  been  made  to  improve  upon  the  previous 
one-point  instrument.  Using  a  specially  designed  collection-lens  system,  off-axis  optical 
aberrations  have  been  minimized  effectively  in  the  probe  volume.  Furthermore,  additional 
electronics  have  been  installed  which  permit  absolute  species  concentration  measurements.  As  a 
preliminary  study,  two-point  OH  time-series  measurements  were  taken  in  three  standard 
turbulent  non-premixed  flames  (H3,  DLR  A  and  DLR_B)  [12].  Two-point  statistics,  including 
space-time  correlations,  spatial  autocorrelation  functions  and  integral  length  scales  were 
computed  to  elucidate  the  flame  structure  and  turbulence-chemistry  interactions. 


Understanding  complex  physical  and  chemical  processes  in  turbulent  combustion 
requires  temporally  and  spatially  resolved  measurement  techniques.  Non-premixed  jet  flames 
have  been  studied  extensively  because  of  their  relative  simplicity  and  well-defined  flow  field. 
Measurements  of  minor  species  concentrations  in  turbulent  flames  are  important  owing  to  their 
significance  to  chemical  reactions  and  pollutant  formation.  Hydroxyl,  in  particular,  is  a  good 
flame  marker,  whose  fluctuations  reveal 
important  interactions  between  fluid  mixing 
and  chemical  reactions. 

Time-series  measurements  of 
important  scalars,  such  as  temperature  [13] 
and  concentrations  of  OH  [6]  and  CH  [5], 
have  been  demonstrated  using  high 
repetition-rate,  laser-based  techniques. 

These  measurements  provide  spectral 
characteristics  and  relevant  time  scales  for 
scalars;  however,  as  pointwise 
measurements,  they  give  no  information  on 
spatial  structures.  To  circumvent  this 
limitation,  simultaneous  multipoint  or  even 
ID  or  2D  time-series  measurements  are 
required.  Examples  of  this  endeavor  include 
multipoint  OH  chemiluminescence  measurements  [14],  two-point  temperature  measurements 
[15],  and  high  repetition-rate  planar  laser-induced  fluorescence  (PLIF)  [16].  For  minor  species, 
two-point  PITLIF  was  previously  developed,  which  provided  quantitative  concentration 
measurements  of  radical  species  at  relevant  turbulent  temporal  and  spatial  scales  [10,  11]. 

Modifications  have  been  made  to  improve  the  earlier  version  of  the  two-point  PITLIF 
system.  Optical  aberrations  were  found  to  have  significant  effects  on  our  previous  two-point 
statistics  [9].  In  particular,  aberrations  occurred  when  either  probe  volume  deviated  from  the 
optical  axis  of  the  detection  system.  Furthermore,  at  large  displacement  distances,  partial 
overlap  arose  between  images  of  the  two  probe  volumes.  This  overlap  created  artificial 


Displacement  (mm) 

Fig.  4.  Spatial  autocorrelation  function  of  H3 
flame  taken  at  the  location  of  peak  [OH]  at 
x/D=20,  using  both  the  old  and  new  two-point 
optical  systems. _ 
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coherence  between  signals  from  the  two  probe  volumes,  which  was  exhibited  as  a  secondary 
peak  in  the  spatial  autocorrelation  function,  as  depicted  in  Fig.  4. 

To  address  this  problem,  we  designed  a  new  lens  system  via  ZEMAX  software,  which 
now  produces  an  aberration-limited  blur  spot  of  less  than  130  pm  at  object  positions  up  to  6  mm 
away  from  the  optical  axis.  As  shown  in  Fig.  5,  the  current  optical  system  is  composed  of  one 
custom-made  aspherical  lens  (LI),  a  10-cm  diameter  UV  beamsplitter  (BS),  and  two  pairs  of 
plano-convex  and  meniscus  lenses  (L2  and  L3;  L4  and  L5).  The  new  setup  has  a  magnification 
ratio  of  1.19  and  a  working  fW  of  2.8.  Each  probe  volume  is  485x485x210  pm\  which  is 
determined  by  the  diameter  of  the  laser  beam  and  the  entrance  slit  width  of  the  monochromators. 
P1TL1F  utilizes  the  principle  of  triple-bin  integration  [10]  which  can  correct  for  the  effect  of 
electronic  quenching  on-the-fly.  The  insertion  of  four  additional  multichannel  scalers  allows  for 
absolute  species  concentration  measurements  at  both  probe  volumes. 

To  demonstrate  the  effectiveness  of  this  revised  setup,  two-point  PITLIF  data  were  taken 
at  the  same  position  (peak  [OH]  at  x/D=20)  within  the  same  flame  used  by  Zhang  et  al  [9].  The 
derived  spatial  autocorrelation  function  shown  in  Fig.4  displays  an  approximately  exponential 
decay  without  a  secondary  peak,  as  expected  from  the  typical  statistics  of  a  turbulent  field. 

For  this  study,  three  non-premixed  flames  (H3,  DLR  A  and  DLR  B)  were  investigated 
chosen  from  the  Turbulent  Nonpremixed 
Flames  (TNF)  Workshop.  Two-point  OH 
measurements  were  taken  at  locations 
corresponding  to  the  mean  peak  [OH]  as 
well  as  positions  along  the  jet  centerline, 
with  x/D  ranging  from  10.00  to  53.25, 
where  x/D=35  and  x/D=53.25  correspond 
to  the  respective  flame  tips  of  the  H3  and 
DLR  A  flames.  The  mean  OH 
concentrations  were  compared  with  the 
available  data  from  the  TNF  database 
[12]  (the  detailed  comparison  can  be 
found  in  Zhang  et  al  [10]).  The 
measurements  agree  quite  closely  except 
at  two  locations:  the  jet  centerline  at 
x/D=50  in  DLR  A  and  peak  [OH]  at 
x/D=20  in  DLR  B.  This  discrepancy 
likely  is  due  to  deviations  in  position 
between  the  two  measurements  and 
increased  uncertainty  at  the  jet  centerline  owing  to  the  small  signal-to-background  ratio. 

Given  two  simultaneous  [OH]  time  series,  a  number  of  two-point  statistics  can  be 
evaluated,  including  the  space-time  correlation,  spatial  autocorrelation  function,  and  integral 
length  scale.  The  normalized  space-time  correlation  is  defined  as  [1 1], 

(OH  (r,t)i)H  (r  +  A r,t+  At)) 

(oh  (r,t  j(oH  (r  +  A r,t  +  A 


f st  (Ar ,  A/) 


LI  BS  L2  L3 


Fig.  5.  Schematic  of  detection  optics  in  the  current  two- 
point  PITLIF  system.  L,  lens;  BS,  beam  splitter;  DET, 
detector  assembly  consisting  of  a  monochromator  and 
photomultiplier  tube.  Solid  and  dashed  rays  identify  the 
signal  from  two  different  probe  volumes. _ 
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where  OH'{r,t)  is  the  fluctuating  component  of  hydroxyl  concentration  at  radial  location  r  and 
time  t,  while  A r  and  At  are  the  displacement  between  the  two  probe  volumes  and  the  lag  time, 
respectively.  Figure  6  displays  the  calculated  space-time  correlation  at  two  peak  [OH]  locations 
in  the  H3  flame.  Positions  near  the  flame  tip  show  a  broader  peak  along  the  displacement  axis, 
indicating  a  larger  coherence  region. 


-u  z.  ■  1  ~ i - i - ' - i- -1 — ■ — — i 

-8-6-4  -2  02  46  8 

Displacement  (mm) 

Fig.  7.  Normalized  spatial  autocorrelation 
functions  for  the  DLR  A  flame:  (a)  at  peak  [OH] 
locations. 


Fig.  6.  Space-time  correlation  at  two  peak  [OH] 
locations  in  the  H3  flame:  (a)  x/D=10  (shaded 
surface);  (b)  x/D=30  (mesh). _ 


The  normalized  spatial  autocorrelation  function  is  related  to  the  space-time  correlation  by 


/.(*••  A'=°> 


(10) 


Figure  7  displays  spatial  autocorrelation  functions  for  the  DLR  A  flame  at  peak  [OH]  locations. 
The  full-width-at-half-maximum  (FWHM)  of  the  curve  becomes  greater  farther  downstream, 
which  is  a  manifestation  of  the  broadening  coherence  region  mentioned  earlier.  In  general,  the 
spatial  autocorrelation  function  is  symmetric  between  the  air  and  fuel  sides.  Exceptions  occur 
primarily  near  the  nozzle  (x/D  <  10).  This  abnormality  might  be  explained  by  the  evolution  of 


Fig.  8.  Hydroxyl  integral  length  scales  for  H3 
flame. 
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Fig.  9.  Hydroxyl  integral  length  scales  for  DLR  A 
and  DLR  B  flames. 


flame  structures.  In  turbulent  jet  diffusion  flames,  the  reaction  zone  lower  in  the  flame  (e.g. 
x/D<5  for  the  DLR  A  flame)  is  composed  of  thin  layers  located  outside  the  turbulent  fuel  jet 
[18].  Because  of  relaminarization,  the  reaction  zone  essentially  acts  as  a  shield  blocking  the 
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penetration  of  fuel  vortices;  hence,  the  reaction  zone  and  air  side  remain  laminar.  Farther 
downstream,  the  flame  front  becomes  more  wrinkled  as  the  effect  of  laminarization  diminishes; 
therefore,  the  asymmetry  of  the  spatial  autocorrelation  function  at  x/D=10  probably  is  due  to 
residual  laminarization. 

The  OH  integral  length  scale  can  be  determined  by  integrating  the  spatial  autocorrelation 
function  over  its  displacement,  A r,  so  that 

(11) 

While  the  spatial  autocorrelation  function  exhibits  asymmetry  near  the  nozzle,  the  calculated 
integral  length  scales  are  about  the  same  on  both  sides  (Fig.  8).  The  hydroxyl  length  scales  are 
well  resolved  and  increase  linearly  with  axial  distance  at  peak  [OH]  locations  (Figs.  8  and  9), 
consistent  with  observations  of  integral  length  scales  for  mixture  fraction  at  mean  stoichiometric 
radii  in  methane/air  jet  diffusion  flames  [19].  This  result  suggests  that  OH  fluctuations  at  peak 
[OH]  locations  essentially  are  dominated  by  flame  motions,  a  confirmation  of  flamelet  theory. 
The  integral  length  scales  basically  are  constant  along  the  jet  centerline  below  the  flame  tip  and 
increase  rapidly  above  the  flame  tip. 


Axial  Measurement  System  and  Results 

The  previously  reported  two-point  OH  time  series  for  turbulent  jet  flames  were  obtained 
with  two-point  separation  along  the  radial  direction  of  the  jet  (Ar),  i.e.,  the  propagation  direction 
of  the  horizontal  laser  beam.  For  a  thorough  understanding  of  OH  fluctuations,  we  also  desire 
two-point  statistics  using  separations  along  other  directions.  As  the  jet  configuration  is 
statistically  axisymmetric  so  that  statistical  analyses  along  the  azimuthal  direction  can  be 
considered  insignificant,  the  only  remaining  direction  of  interest  is  along  the  axis  of  the  jet,  i.e., 
the  z-direction.  Investigation  of  axial  two-point  statistics  requires  a  different  laser  configuration 
than  that  of  the  previous  radial  two-point  system:  either  a  single  vertical  laser  beam  or  two 
horizontal  beams  with  adjustable  vertical  separation.  The  combustor  used  in  this  study  has  a 
closed  base  that  does  not  permit  vertical  alignment  of  a  laser;  hence,  the  two-beam  configuration 
has  been  adopted.  With  the  installation  of  our  new  diode-pumped  laser  system,  which  generates 
nearly  ten  times  more  UV  power  than  before,  the  laser  now  can  be  split  into  two  beams  while 
still  providing  a  sufficient  signal  level  at  each  probe  volume.  This  section  describes 
experimental  details  for  this  new  approach  and  its  demonstration  within  a  well-studied  turbulent 
hydrogen-nitrogen  nonpremixed  flame,  the  H3  flame  (H2/N2  =  50:50,  Re  =  10,000). 

For  this  variation  on  the  PITLIF  technique,  the  mode-lock  Ti:  Sapphire  laser  is  pumped  by 
a  532-nm  diode-pumped  cw  laser  at  13.5  W.  The  generated  IR  pulsed  beam  at  a  repetition  rate 
of  80  MHz  then  is  frequency-tripled  into  the  UV.  The  output  UV  laser  pulse  has  a  temporal 
width  of  2  ps  and  a  spectral  width  of  0.3  nm.  Its  wavelength  is  tuned  to  306.5  nm  to  excite 
multiple  transitions  in  the  (v'  =  0,  v"  =  0)  vibrational  band  of  the  A  2S+  -  X  2IIj  electronic 
system  of  OH  selectively.  The  excitation  wavelength  is  selected  to  minimize  the  temperature 
dependence  of  the  overall  Boltzmann  fraction,  so  that  quenching-corrected  fluorescence  signals 
can  be  converted  directly  to  OH  concentrations.  Simulations  of  the  OH  absorption  spectrum 
versus  temperature,  using  the  spectroscopic  code  developed  by  Seitzman  [20],  verify  a  variation 
in  absorption  coefficient  of  only  ±5%  at  1500-2500  K. 
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Figure  10.  Schematic  of  the  laser  alignment  used  for  the  current  two-point  PITLIF  system: '/« X,  zero-order 
quarter  waveplate;  M,  mirror;  BS,  polarizing  cube  beamsplitter;  BC,  polarizing  cube  beam  combiner;  PV,  probe 
volume;  L,  lens;  BD,  beam  dump.  The  red  and  green  arrows  identify  the  polarization.  “P”  and  “S”  designate 
two  split  laser  beams  based  on  their  polarizations. 


Figure  10  shows  the  experimental  layout  used  for  aligning  the  two  beams  of  this 
investigation.  The  linearly  polarized  UV  beam  from  the  laser  system  first  is  converted  into  a 
circularly  polarized  beam  with  a  zero-order  quarter  waveplate  (V*  X.)  and  then  split  into  two 
beams  (“P”  and  “S”)  with  different  polarizations  within  a  polarizing  beamsplitter  (BS).  The  two 
split  beams  subsequently  are  folded  with  mirrors  (M2  and  M3)  and  recombined  with  a  beam 
combiner  (BC)  to  form  two  parallel  horizontal  beams  across  the  flame.  The  axial  two-point 
separation  then  can  be  adjusted  continuously  by  translating  the  M3  mirror.  To  ensure  day-to-day 
consistency  in  performance,  all  these  optical  components  except  the  wave-plate  and  lens  are 
mounted  on  a  30  cm  x  30  cm  vertically  placed  optical  breadboard.  This  arrangement  of  two 
horizontal  laser  beams  offers  much  freedom  for  two-point  measurements,  as  the  separation 
between  the  two  probe  volumes  (PV)  can  be  made  radial,  axial,  or  both. 

The  optical  detection  system  and  the  subsequent  photon-counting  system  are  essentially 
the  same  as  those  used  in  previous  radial  two-point  experiments  [21].  To  summarize, 
fluorescence  photons  from  two  probe  volumes  are  collected  by  an  aspherical  lens,  split  by  a 
beamsplitter,  and  eventually  focused  by  lens  doublets  onto  corresponding  detector  assemblies, 
with  each  consisting  of  a  monochromator  and  PMT.  Adjustment  of  two-point  separation  is 
achieved  by  translating  one  detector  while  keeping  the  other  one  fixed.  As  the  entrance  slit  of 
each  monochromator  defines  the  location  of  the  corresponding  probe  volume,  calibration  is 
required  to  establish  the  ratio  between  translation  of  the  probe  volume  and  that  of  the  detector  by 
imaging  a  needle  tip  placed  at  the  probe  volume.  For  both  the  radial  and  axial  directions,  this 
ratio  is  found  to  be  1.19,  which  matches  the  magnification  ratio  of  the  collection  lens  system. 
Each  probe  volume  has  dimensions  of  350x350x210  pm3,  as  defined  by  the  laser-beam  diameter 
(e'2)  and  the  entrance-slit  width  of  the  monochromator. 

A  new  problem  that  emerges  for  axial  two-point  measurements  is  interference  between 
the  two  laser  beams.  The  p-  and  s-  polarized  beams  have  a  vertical  dimension  of  330  and  190 
pm  at  the  focal  point,  respectively,  which  implies  full  or  partial  overlap  between  the  two  beams 
for  A z  ranging  between  -260  and  260  pm;  moreover,  the  height  of  the  entrance  slits  for  the 
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monochromators  was  set  at  approximately  750  pm,  corresponding  via  magnification  effects  to 
625  pm  of  vertical  acceptance  window  at  the  probe  volumes.  This  fact  suggests  that  it  is  still 
possible  to  collect  fluorescence  photons  from  both  laser  beams,  even  though  the  two  beams  are 
well  separated.  To  avoid  all  interference  from  the  two  beams,  the  minimum  two-point  separation 
|Az|  should  be  around  475  pm.  Experimentally,  the  two-beam  interfering  range  for  A z  can  be 
determined  by  examining  the  change  in  signal  detected  by  the  fixed  PMT  versus  two-point 
separation.  This  interference  range  can  be  defined  as  that  beyond  which  the  photon  rate  has 
minimal  variation.  For  all  measurements,  the  interference  range  was  found  to  be  between  -1.0 
and  0.5  mm;  hence,  measurements  within  this  A z  range  were  rejected.  For  the  systems  studied  so 
far,  this  minimum  separation  has  little  effect  on  our  ability  to  resolve  coherent  structures  within 
turbulent  flames. 

Space-time  correlations  derived  from  the  axial  two-point  time-series  measurements  are 
shown  in  Fig.  1 1  for  z/D  =10.  At  zero  two-point  separation  (Az  =  0),  the  space-time  correlation 
reduces  to  the  temporal  autocorrelation  function,  from  which  we  may  determine  the  integral  time 
scale.  This  measurement  thus  provides  an  independent  way  of  measuring  this  macroscale,  which 
is  found  to  match  within  10%  the  results  from  previous  single-point  and  radial  two-point 
measurements  [22].  Compared  with  the  previously  reported  radial  space-time  correlations,  the 
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Figure  1 1 .  Space-time  correlations  of  [OH]  fluctuations  for  two-point  separations  along  the  axial  direction  at 
the  radial  peak  [OH]  location  corresponding  to  z/D  =  1 0.  Negative  Az  corresponds  to  the  upstream  direction 
and  vice  versa. 


axial  space-time  correlations  have  some  peculiar  features.  First,  for  nonzero  two-point 
separations  at  upstream  locations,  i.e.,  when  the  second  probe  volume  is  moved  upstream,  each 
correlation  curve  has  its  maximum  at  a  certain  time  delay.  This  feature  appears  because  OH 
structures  are  convected  from  upstream  locations  and  take  this  amount  of  time  to  reach  the  fixed 
probe  volume.  Hence,  the  time  delay  is  equal  to  the  two-point  separation  divided  by  the 
convection  velocity.  Second,  a  progressive  decrease  below  unity  occurs  in  the  correlation 
maxima  with  increasing  two-point  separation.  The  reason  for  this  behavior  is  that  during  the 
convection  process,  the  OH  structures  experience  chemical  reactions  and  turbulent  mixing,  so 
that  the  longer  the  time  it  takes  to  reach  the  fixed  probe  volume,  the  less  the  structures  look  like 
the  original  ones.  A  similar  phenomenon  has  been  observed  for  space-time  correlations  of 
velocity  in  a  nonreacting  jet  (Wills)  [23]  and  for  those  of  temperature  in  premixed  jet  flames 
(Ghenai  and  Gokalp)  [24], 
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Convection  of  turbulence  and  scalar  structures  stems  from  Taylor’s  original  hypothesis 
on  grid  turbulence.  According  to  Taylor  [25],  at  low  turbulence  levels  the  spatial  pattern  of 
fluctuating  velocity  is  rigidly  translated  to  other  locations  at  the  mean-flow  velocity;  hence, 


fsr  -  VcTv^)=  fsr  (Az>A/  +  r,),  (12) 

where  Vc  is  the  convection  velocity. 

This  hypothesis,  the  so-called  frozen- 
turbulence  approximation,  is  accurate 
only  for  grid  turbulence  and  typically  is 
not  strictly  valid  for  free  shear  flows, 
such  as  jets  (Pope)  [26].  Space-time 
correlations  measured  in  this 
investigation  confirm  that  gradual 
distortion  occurs  for  OH  structures 
when  being  convected  downstream 
within  the  jet. 

Even  though  Taylor’s  hypothesis 
is  not  strictly  correct,  convection  of 
turbulent  eddies  and  accompanying 
scalar  structures  at  steady  velocity  is  a 
prevalent  concept.  The  convection 
velocity,  which  may  not  be  necessarily 
the  mean-flow  velocity,  then  can  be 
defined  as  the  two-point  separation 
divided  by  the  time  delay  corresponding  to  the  correlation  maximum, 


Figure  12.  Convection  velocity  of  OH  structures  estimated 
from  space-time  correlations  of  various  two-point  separation 
distances  at  the  radial  peak  [OH]  location  for  z/D  =  10.  The 
effective  convection  velocity  is  defined  as  the  average  (the 
dashed  line)  after  discarding  the  first  and  last  three  points 
owing  to  errors  from  truncation  and  curve-fitting. 


where  the  delay  time  zd  can  be  described  mathematically  as 


(13) 

(14) 


A  cubic  spline  curve  can  be  used  to  fit  the  correlation  curve  fsi{Az,  At)  to  find  the  delay 
time,  id,  through  Eq.  (12).  The  calculated  convection  velocity  at  the  radial  peak  [OH]  location 
for  z/D  =  10  is  shown  in  Fig.  12.  Clearly,  the  convection  velocity  is  not  a  single  value,  and 
instead  depends  somewhat  on  the  choice  of  two-point  separation.  This  ambiguity  in  convection 
velocity  has  been  noted  previously  for  velocity  measurements  in  non-reacting  jets  [23];  variable 
choices  were  proposed  for  this  velocity  that,  however,  do  not  shed  more  light  on  the  underlying 
physics.  For  this  investigation,  an  effective  convection  velocity  is  determined  as  the  average  of 
convection  velocities  at  different  two-point  separations.  Points  are  discarded  at  small  and  large 
separations,  as  such  calculations  are  likely  contaminated  owing  to  truncation  errors  at  small 
separations  and  errors  resulting  from  curve-fitting  at  large  separations.  We  note  that  other  ways 
are  available  for  calculating  the  convection  velocity,  either  by  the  constant  time-delay  method 
[23]  or  based  on  cross-spectra  [27],  which  may  offer  greater  significance  for  acoustic  noise 
generation. 
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The  space-time  correlation  for  downstream  two-point  separations  in  Fig.  12  is  quite 
different  from  that  for  upstream  separations.  The  maximum  of  each  correlation  curve  occurs  at 
zero  delay  time,  and  undergoes  faster  decay  with  respect  to  increasing  two-point  separation.  The 
reason  for  this  behavior  is  that  downstream  structures  cannot  travel  back  and  thus  are  unable  to 
relate  better  to  the  fixed  probe  volume  for  At  >  0.  This  feature  illustrates  that  the  space-time 
correlation  for  a  stationary  process  is  typically  not  symmetric  with  spatial  separation;  rather, 

fsr  (Az,A/)=  /„.  (-Az,-A/)*  (-Az,  A/),  (15) 

which  can  be  shown  easily  by  switching  the  roles  of  the  two  probe  volumes.  The  different 
features  of  axial  space-time  correlations 
between  positive  and  negative  separations  also 
should  be  observed  for  radial  space-time 
correlations  at  non-zero  radial  positions; 
however,  owing  to  the  much  smaller  radial 
velocity  for  axisymmetric  jets,  the  radial 
convection  of  scalar  structures  is  much  less 
prominent,  which  may  explain  the  failure  to 
observe  such  convection  for  the  thermal 
structures  of  (Ghenai  and  Gokalp)  [24]  . 

A  different  way  of  studying  the  space- 
time  correlation  is  to  develop  a  constant 
correlation  map.  For  an  ideal  Taylor’s  model, 
constant  space-time  correlations  satisfying  Eq. 

(9)  would  be  shown  in  the  (Az,  A/)-plane  as 
parallel  lines  with  a  slope  equal  to  the 
convection  velocity,  -Vc,  indicating  rigid 
convection  of  scalar  structures  at  a  constant 
mean  velocity.  As  shown  in  Fig.  13,  space- 
time  correlation  maps  of  OH  concentration 
retain  some  features  of  Taylor’s  model.  The 
constant  correlation  lines  are  approximately 
parallel  to  each  other  over  a  narrow  range  of 
two-point  separations;  however,  they 
eventually  form  close  loops  over  the  large, 
negative  Az  region,  further  confirming  the  idea 
of  non-rigid  convection  of  scalar  structures 
within  jets.  Furthermore,  the  variable  slopes  of 
the  constant  correlation  lines  make  the 
definition  of  an  unambiguous  convection 
velocity  impossible;  instead,  the  convective 
velocity  becomes  a  function  of  spatial 
separation,  as  indicated  by  Fig.  13.  Compared 
with  the  radial  peak  [OH]  location  at  z/D  =  10, 
the  space-time  correlations  at  the  jet  axis  for 
z/D  =  45  offer  greater  resemblance  to  an  ideal 
Taylor’s  model.  A  survey  over  all  measurements  of  this  study  further  suggests  that  space-time 
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Figure  13.  Space-time  correlation  maps  of  OH 
concentration  fluctuations  for  two-point 
separations  along  the  axial  direction  at  radial 
peak  [OH]  for  z/D  =  10  and  at  the  jet  axis  for  z/D 
=  45.  Negative  Az  corresponds  to  the  upstream 
direction. 
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correlations  at  the  jet  axis  and  at  higher  z/D  more  closely  approach  Taylor’s  ideal  case  than  those 
off  the  centerline  and  at  lower  z/D.  This  phenomenon  is  quite  strange,  as  axial  positions  at  high 
z/D  have  higher  turbulence  levels  so  that  Taylor’s  hypothesis  should  be  more  likely  to  fail. 

Conclusions  and  Recommendations 

A  two-point  high-speed  laser-induced  fluorescence  technique  (PITL1F)  for 
simultaneously  capturing  fluctuations  of  minor  species  at  two  spatial  locations  has  been 
developed.  Various  issues  in  instrumentation  for  two-point  PITLIF  have  been  presented  and 
applications  of  this  technique  have  been  demonstrated  for  measurements  of  OH  in  a  variety  of 
turbulent  nonpremixed  and  partially  premixed  jet  flames. 

Conclusions 

The  following  conclusions  can  be  made  from  this  investigation. 

(1)  Simultaneous  two-point  quenching-corrected  measurements  have  been  made  for  the  first 
time  at  sampling  rates  reaching  40  kHz.  This  system  offers  sufficient  spatial  and 
temporal  resolution  for  investigating  two-point  statistics  of  minor  species  concentrations, 
including  the  space-time  correlation,  spatial  autocorrelation  and  integral  length  scale.  In 
principle,  this  technique  can  be  applied  to  a  variety  of  minor  species,  including  OH,  CH, 
and  NO,  provided  that  enough  laser  power  can  be  obtained  at  the  excitation  wavelength 
of  the  molecule  of  interest. 

(2)  Optical  aberrations  of  the  optical  detection  system  are  important  for  both  single-point  and 
two-point  measurement  techniques.  Such  aberrations  have  a  significant  impact  on  spatial 
resolution  and  the  measured  two-point  statistics;  hence,  meticulous  optical  design  should 
be  performed  to  minimize  optical  aberrations. 

(3)  For  all  flames  of  this  investigation,  fluctuations  in  fluorescence  lifetime  have  a  negligible 
effect  on  the  derived  two-point  statistics.  Hence,  two-point  statistics  calculated  from 
fluorescence  bin  counts  can  be  treated  as  concentration;  therefore,  uncertainties  involved 
in  the  quenching  correction  process  can  be  avoided  and  high  sampling  rates  can  be 
achieved,  which  provides  the  potential  for  studying  flames  at  high  Reynolds  numbers. 

(4)  For  turbulent  nonpremixed  jet  flames,  the  temporal  autocorrelation  functions  exhibit 
universal  self-similarity,  while  the  radially  spatial  autocorrelation  functions  are  self¬ 
similar  for  small  to  medium  two-point  separations,  owing  to  the  non-universality  of 
large-scale  fluctuations  that  suggests  the  hydroxyl  integral  time  scale  is  sufficient  for 
characterizing  macroscale  temporal  fluctuations;  however,  the  integral  length  scale  may 
be  meaningful  only  to  a  certain  spatial  extent,  for  instance,  within  the  corrugation  scale  of 
the  OH  structures.  The  existence  of  coherent  vortices  near  the  flame  base  is  found  to  be 
responsible  for  negative  spatial  correlations,  which  become  less  appreciable  with 
increasing  jet  Reynolds  number. 

(5)  For  turbulent  nonpremixed  jet  flames,  the  radial  hydroxyl  integral  length  scale  displays 
different  trends  at  mean  radial  [OH]  peaks  and  close  to  the  flame  tip.  At  mean  radial 
[OH]  peaks,  the  OH  integral  length  scale  is  related  to  corrugation  of  the  OH  layer  and  is 
dictated  by  the  average  eddy  size,  thus  following  a  linear  relation  with  axial  distance. 
Around  the  flame  tip,  the  OH  integral  scale  is  close  to  that  of  a  conserved  scalar  and  of 
the  dynamic  scale.  Below  the  flame  tip,  the  OH  integral  length  varies  little  with 
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increasing  z/D  owing  to  the  effect  of  enhanced  dissipation  of  large-scale  vortices  at 
increasingly  high  temperatures.  Beyond  the  flame  tip,  the  OH  length  scale  increases 
more  than  linearly  with  z/D  because  of  penetration  of  large-scale  buoyancy-induced 
vortices.  The  radial  integral  length  scale  depends  weakly  on  the  jet  Reynolds  number  but 
strongly  on  fuel  dilution  owing  to  its  influence  on  flame  length  and  overall  flame 
structures. 

(6)  Hydroxyl  integral  time  scales  in  turbulent  nonpremixed  jet  flames  can  be  scaled 
reasonably  by  the  ratio  of  OH  integral  length  scale  and  the  primary  velocity  component 
along  the  same  direction.  The  integral  time  scale  depends  strongly  on  both  the  jet 
Reynolds  number  and  fuel  dilution.  For  the  flames  under  investigation,  OH  integral  time 
scales  can  be  related  approximately  linearly  to  the  mean  flame-front  passing  period  at 
most  locations  within  flames,  suggesting  that  the  integral  length  scale  indicates  the 
mobility  of  the  flame  front. 

(7)  For  turbulent  partially  premixed  jet  flames,  temporal  and  radially  spatial  autocorrelation 
functions  are  generally  not  self-similar,  suggesting  that  the  OH  integral  time  and  radial 
length  scales  cannot  be  treated  as  parameters  characterizing  solely  macroscale 
fluctuations  of  [OH]  in  these  flames.  At  sufficient  premixing  levels,  the  flames  exhibit  a 
double-flame  structure,  where  the  premixed  front  displays  more  rapid  [OH]  fluctuations 
than  the  nonpremixed  front.  The  spatial  autocorrelation  function  displays  strong  negative 
values,  including  at  far  downstream  locations,  suggesting  the  existence  of  a  prolonged 
region  of  coherent  vortices  and  laminarization,  which  may  result  from  the  thick  high- 
temperature  region  between  the  two  flame  fronts. 

(8)  For  turbulent  partially  premixed  jet  flames,  temporal  fluctuations  of  [OH]  at  the  premixed 
flame  front  are  characterized  by  both  a  fast  time  scale  corresponding  to  fast  flamelet 
crossings  and  a  slow  time  scale,  which  is  approximately  ten  times  smaller  than  the  fast 
scale.  At  the  nonpremixed  front,  the  integral  time  scale  follows  a  scaling  relation  similar 
to  that  for  purely  nonpremixed  jet  flames. 

(9)  For  turbulent  partially  premixed  jet  flame,  the  radial  integral  length  scale  at  the  premixed 
flame  front  is  greater  than  that  at  the  nonpremixed  flame  front  for  the  same  axial  distance. 
At  the  nonpremixed  front,  the  radial  integral  length  scale  varies  in  a  nonlinear  fashion 
with  respect  to  axial  distance,  which  may  be  related  to  the  prolonged  laminarized  region. 

(10)  Axial  two-point  OH  statistics  for  nonpremixed  flame  exhibit  prominent  features  of  OH 
convection.  Convection  of  OH  structures  occurs  at  approximately  the  flow  velocity  and 
deviates  from  Taylor’s  hypothesis  with  distortion  of  OH  structures.  The  axial  length 
scale  is  generally  greater  than  the  radial  length  scale,  again  suggesting  elongated  OH 
structures,  which  become  isotropic  at  far  downstream  locations.  This  integral  length 
scale  does  not  follow  the  dynamic  length  scale,  as  measured  in  turbulent  nonreacting  jet, 
which  suggests  the  influence  of  heat  release  and  localization  of  OH  structures. 

In  general,  OH  chemistry  does  not  seem  to  have  an  appreciable  effect  on  macroscale 
statistics,  including  space-time  correlations  and  integral  time  and  length  scales.;  however, 
such  chemistry  may  contribute  to  high-frequency  components  of  OH  PSDs,  leading  to  the 
absence  of  a  power-law  region,  as  would  be  expected  for  a  scalar  in  a  turbulent  system. 

A  method  for  studying  the  statistics  of  flame  front  motion  has  been  devised,  by  threshold- 
filtering  the  OH  concentration.  A  comparison  between  flame-front  and  OH  statistics 
shows  that  these  two  statistics  are  comparable  over  most  measurement  locations. 
However,  deviations  occur  near  the  flame  tip,  where  OH  structures  are  typically  thick. 
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An  ambiguous  definition  of  the  flame  front  in  this  region  is  challenging,  as  the 
combustion  zone  is  typically  thin  throughout  a  flame  [18].  The  validity  of  this  analysis  is 
subjected  to  further  scrutiny. 

The  distinction  between  partially  premixed  and  nonpremixed  jet  flames  is  interesting.  A 
prolonged  laminarized  region  resulting  from  the  high  temperature  between  the  two  flame 
fronts  is  proposed  as  responsible  for  the  peculiar  spatio-temporal  statistics  of  partially 
premixed  flames.  Possible  crossing  of  the  premixed  front  at  the  mean  nonpremixed  peak 
may  result  in  the  measured  bimodal  PDF  in  OH  concentration.  Other  issues  also  require 
further  clarification,  such  as  the  slow  time  scale  of  OH  fluctuations  at  the  premixed  front 
and  the  cause  of  the  visible  flame-tip  pulsation. 

Recommendations 

Based  on  the  progress  to  date,  the  following  recommendations  can  be  given. 

(1)  The  current  two-point  PITLIF  setup  can  be  improved  further  by  installing  separate 
detection  optics  for  the  two  different  detectors,  with  the  detection  optics  and  detector 
corresponding  to  the  movable  probe  volume  mounted  on  a  single  translation  stage.  This 
approach  offers  two  advantages.  First,  it  avoids  off-axis  optical  aberrations,  as  the 
movable  probe  volume  is  always  on  the  optical  axis.  Second,  it  avoids  splitting  of  the 
fluorescence  signal  with  a  beamsplitter,  thus  doubling  the  detected  signal  level.  In 
addition,  the  laser  system,  including  the  Tsunami  laser  and  the  harmonic  generator,  may 
be  put  in  an  enclosure  with  local  climate  control.  The  ensuing  advantage  is  avoidance  of 
dramatic  fluctuations  in  temperature  and  humidity,  which  are  typical  when  running  a 
flame.  This  tactic  is  quite  critical  for  ensuring  stable  laser  operation. 

(2)  Though  a  powerful  technique,  two-point  PITLIF  is  a  single-scalar  measurement,  which 
cannot  reveal  a  full  picture  of  turbulent  flames.  As  processes  in  turbulent  combustion  are 
highly  correlated  in  a  nonlinear  fashion,  it  is  quite  desirable  to  make  multi-scalar 
measurements.  With  the  current  facility,  it  is  possible  to  perform  simultaneous  OH  and 
temperature  measurements,  as  temperature  can  be  measured  by  an  available  argon-ion 
laser.  With  simultaneous  [OH]  and  T,  the  effects  of  scalar  dissipation  rate  can  be 
investigated  for  both  scalars  in  turbulent  nonpremixed  flames,  which  is  valuable  for 
turbulent  combustion  modeling.  As  T  can  be  related  solely  to  mixture  fraction  in  a  fuel- 
lean  or  fuel  rich  condition,  we  also  may  use  two-point  [OH]  measurements  to  define  the 
stoichiometric  condition  by  tracking  the  movement  of  the  flame  front.  The  simultaneous 
T  time  series  may  be  converted  to  a  mixture  fraction  time  series,  which  provides  more 
significance  than  a  species  time  series.  Additionally,  the  simultaneous  [OH]  and  T  data 
definitely  will  shed  more  light  on  less  understood  turbulent  partially  premixed  flames. 

(3)  The  current  technique  may  be  extended  to  applications  involving  CH  or  NO.  For  NO,  the 
quadrupled  output  of  the  harmonic  generator  can  be  used  (-226  nm).  By  doping  the  flow 
with  NO,  NO  concentration  can  be  used  as  a  surrogate  for  studying  mixing  processes 
within  turbulent  nonreacting  flows  or  even  flames  for  which  NO  production  and 
destruction  are  less  significant  as  compared  to  any  dopant  concentration.  For  CH,  the 
doubled  output  of  the  harmonic  generator  can  be  used  (-431  nm).  As  CH  is  a  better 
marker  of  the  flame  front  than  OH  [18],  its  time  series  may  provide  better  pictures  of 
turbulence-flame  interactions. 

(4)  Time  and  space  series  (TASS)  simulations  can  be  adopted  to  extend  the  use  of  two-point 
OH  time-series  measurements  in  turbulent  nonpremixed  jet  flames.  With  this  method, 
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realistic  mixture  fraction  time  series  can  be  simulated,  along  with  integral  time  and  length 
scales  of  mixture  fraction. 

Finally,  with  two-point  PITLIF  producing  rich  information  on  OH  fluctuations,  it  is 
desirable  to  collaborate  with  combustion  modelers  to  make  use  of  these  data  to  validate 
advanced  combustion  models. 
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